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Abstract 

At present there is good agreement between the neutrino mass-squared difference determined 
from the solar neutrino data and the anti-neutrino mass-squared difference determined from the 
KamLAND reactor anti-neutrino experiment. However, the central values of the two cases differ 
from each other by about 10 -5 eV 2 . An improvement in the accuracy of both the solar neu- 
trino experiments and reactor anti-neutrino experiments can establish the existence of a non-zero 
difference between neutrino and anti-neutrino mass-squared differences and provide a signal for 
CPT violation. In this paper, we show how such a difference can arise through the CPT violating 
neutrino mass terms from Planck scale physics. 
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I. INTRODUCTION 



KamLAND experiment recently detected the distortion due to oscillations in the anti- 
neutrino spectrum from the reactors and determined the corresponding mass-square differ- 
ence, A 2 i, to great precision jlj|. At present there is good agreement between A 2 i and the 
mass-squared difference of the neutrinos, A21, determined from the analysis of solar neu- 
trino data P). However, the best-fit values of the two As differ from each other by 10~ 5 
eV 2 . Together, solar and KamLAND data impose the constraint |A 2 i — A 2 i| < 1.1 x 10~ 4 
eV 2 [3]. Future reactor experiments, located at a distance of about 70 Km from the source 
so that the oscillation minimum coincides with spectral maximum, are expected to improve 



the precision of A 21 even further 
LENS P| and other experiments |6| 



Similarly future solar neutrino experiments, such as 
are expected to improve the accuracy of A 2i . These 
future experiments may indeed show that there is a non-zero difference between A 21 and 
A 2 i, thus establishing a signal for CPT violation in the neutrino sector 0,0]. 

If A 21 and A 21 are indeed found to be different, a natural question to ask is: How does 
this CPT violation arise? In this letter, we assume that CPT violation in neutrino sector 
arises due to Planck scale effects. We parametrize these effects in terms of Planck scale CPT 
violating neutrino mass terms and calculate the difference between A 2 i and A 21 arising due 
to these terms. 

We assume that neutrino masses mainly arise due to Grand Unified Theory (GUT) dy- 
namics via see-saw mechanism [2j and these masses are CPT conserving. We further assume 
that the CPT violation arises only at Planck scale and parametrize it by the effective neu- 
trino mass term 

McPT = (1) 
where a and (3 are flavour indices. The mass term for anti-neutrinos will have the opposite 
sign 10]. Since these are effective masses arising from Planck scale effects, they are sup- 
pressed by 1/Mpi. Since these are assumed to be the masses at the low energy scale, the 
electroweak vacuum expectation value (VEV) v = 174 GeV is used to make these terms have 
dimension of mass. In eq. the term A is a 3 x 3 matrix in flavour space whose elements 
are of order 1. We assume that the Planck scale interaction is "flavour blind", i.e. the 
elements of the matrix A are independent of a,f3 indices. In this case, the contribution to 
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the neutrino mass matrix is of the form: 



1 1 1 

V 1 1 1 / 



(2) 



where the scale fi is 



fi = v 2 /2M Pl ~ 10~ 6 eV. (3) 

In our calculations, we take eq. @ as a perturbation to the main part of the neutrino mass 
matrix, that is generated by GUT dynamics. 



II. CALCULATION 

The theoretical framework, in which the Planck scale mass terms are treated as pertur- 
bation to GUT scale neutrino masses, is developed in ref. Here we briefly recapitulate 
some salient features of this framework. In the discussion below the labels a and (3 refer to 
flavour eigenstates and the labels i and j refer to mass eigenstates. The unperturbed (0 th - 
order) neutrino mass matrix A4 is diagonalized by a unitary matrix U to yield the matrix 
M, whose eigenvalues are real and non-negative. As stated before, M. is generated by grand 
unified dynamics and is related to U and M through the relation 



/ u el u e2 u e3 \ 



M = U* M U\ where U 



Ufii U^ 

\ U T i U t2 U t3 , J 

We adopt the usual parameterization 1^|: U e2 jU e \ = tan#i 2 , U^ 3 /U r3 = tan^ 23 and \U e3 
sin 6\ 3 . We include all possible phases in the definition of the neutrino mixing matrix, 

U = diag(e 4/l ) IT R(6 13 ) U* R{6 12 ) diag(e 4 



(4) 



(5) 



The phase 5 appearing in II = diag(e i,5/2 ,l,e- i5/2 ) is the one that enters oscillation probabil- 
ities and leads to CP violation in neutrino oscillations. a« are the so called Majorana phases 
and the matrix with this phases has the form diag(e mi , e M2 , 1). The diagonal phase matrix 
on the left is given by diag(e i ^ 1 , e*^ 2 , e*^ 3 ), which are usually absorbed into the definition of 
the respective charged lepton field. It is possible to rotate away the phases fi, if the mass 
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matrix (|3J) is the complete mass matrix. However, since we are going to add another contri- 
bution to this mass matrix, the phases of the zeroth order mass matrix have an impact 
on the complete mass matrix and thus must be retained. By the same token, the Majorana 
phases which are usually redundant for oscillations have a dynamical role to play now. 

Planck scale effects will add other contributions to the mass matrix. The additional term 
has the form given in eq. (J2J) and with its inclusion the neutrino mass matrix becomes 

M — > M' = M + fi X, (6) 

It is possible that the CPT violating mass terms may be arise not at Planck scale but at 
some scale Mx below Mpi but well above the GUT scale (by an order of magnitude or so). 
In such the perturbation parameter [i = v 2 /2Mx, rather than v 2 /2M p i. 

We now define the hermitian matrix A4''A4' and find its eigenvalues and eigenvectors. 
The differences of pairs eigenvalues give us the modified mass-squared differences and the 
eigenvectors give us the modified neutrino mixing matrix. To first order in the small param- 
eter fj,, we have 

M n M' = M^M +M*fj,\. (7) 

This matrix is diagonalized by a new unitary mixing matrix U' . We denote this diagonal 
matrix, correct to first order in /i, to be M' 2 . Using eq. (J3J), we can rewrite Ai in the above 
equation in terms of the diagonal matrix M. Converting A4' also into its diagonal form M', 
we can rewrite the above equation as 

U'M' 2 U' ] = U(M 2 + m ] M + Mm)U ] with m = \i U* X U. (8) 

From the above equation, it follows that the mixing matrix U', correct to first order in /i, is 
related to zeroth order mixing matrix U by 

U' = U (l + i8d), (9) 

where 56 is a hermitian matrix proportional to \l. 

Substituting the expression for U' from eq. © in eq. ^ we obtain 

M 2 + m f M + Mm = M' 2 + [iS9, M' 2 ] . (10) 

Because M' 2 is diagonal, the diagonal terms of [i89, M' 2 ] in the above equation are zero. 
Therefore, to first order in /x, the mass squared differences AM| = M 2 — M 2 get modified 
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as: 

AM| = AM?- + 2 (MiReM - M.-Refm^]), (11) 

where there is no summation over the repeated indices. The above equation shows the 
correction for neutrino mass-squared difference. Because the Planck scale corrections are 
assumed to be CPT violating, the correction for anti-neutrino mass-squared difference will 
have the opposite sign. 

III. RESULTS 

Note from eq.(|ll|) that the correction term depends crucially on the type of neutrino 
mass spectrum. For a hierarchial or inverse hierachial spectrum the correction is negligible. 
Hence we consider a degenerate neutrino spectrum and take the common neutrino mass to 
be 2 eV, which is the upper limit from the tritium beta decay experiment 

From the definition of the matrix m in eq. (jSJ), we find 

mil = /ie* 2ai (U el e th + + U Tl e lf3 f 

m 22 = /ie i2a2 (U e2 e ih + U^ 2 e ih + U T2 e if3 f . (12) 

The contribution of the terms in the Planck scale correction, Re(m22) — Re(mn), can be 
additive or subtractive depending on the values of the phases a± and a 2 . Similarly the 
magnitudes of Re(m22) and Re(mn) are functions of the phases /j. Thus we try to find a 
combination of these phases which can give rise to a significant difference between neutrino 
and anti- neutrino mass-squared difference. In our calculations, we used 9i 2 = 34°, # 13 = 10°, 
$23 = 45° and 5 cp = 0. We define the percentage correction P to be ratio of the difference 
between the corrected neutrino and anti-neutrino mass-squared difference to the unperturbed 
mass squared difference among the first and the second mass eigenstates as; 

P = 100 x (AMfj - ~AM 21 ) / AM 21 (13) 

In Fig.(l) we plot our results as contours of constant P in the a% — a 2 plane, for /j = 0. 
We note that the maximum difference allowed is about ±5% only. This is because the 
magnitudes of m 22 and mn remain relatively small in the limit /$ = 0. In Figs. (2) we plot 
contours of constant P in the f\ — f 2 plane with f 3 = = a\ = a 2 . When one phase is large 
and the other is small, there is a possibility that all the terms in m 2 2 and m u add up to 
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give a large value for their difference and hence a large value for P. We note that a much 
large change, varying from 30% to —20% is possible. In Figs. (3) and (4) we plot contours of 
constant P in f\ — f 3 plane and in f 2 — fz plane. In generating these plots, we set all other 
phases to zero, as in the case of Fig. (2). 

The effect scales directly as the common neutrino mass, so if one takes the WMAP 
[14I ] bound on the common neutrino mass which is 0.23 eV then there is negligible effect. 
However if there are new flavor blind interactions between GUT and Planck scales, then a 
low common mass is compensated by a lower scale of these new interactions and we still 
can get appreciable effects. In fact one can even say that if the gap between A21 and A21 
persists and future data indicate small neutrino mass, this can be an indication of new CPT 
violating flavor blind interactions below Planck scale. 
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FIG. 1: P as a function of the Majorana phases 01 and a>2 



IV. CONCLUSIONS 

Both solar and reactor data are well accounted for by invoking neutrino oscillations. Fit 
to solar data give a large region for the neutrino mass squared difference in the two flavor 
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FIG. 2: P as a function of the phases f\ and / 2 

parameter space. The fit to reactor data however gives a very strongly constrained anti- 
neutrino mass squared difference. The best fits of the two mass squared differences are 
appreciably different from each other. Further improvement in KamLAND systematics and 
future solar neutrino data may remove this discrepancy. However if this mismatch between 
the best fits persists, then CPT breaking in the neutrino sector will be established. We have 
demonstrated that a flavour blind CPT violating neutrino masses from Planck scale physics 
can nicely accomodate this effect. This effect is crucially dependent on the neutrino mass 
spectrum and gives rise to observable difference between A 2 i and A 2 i only for a degenerate 
neutrino mass spectrum with m v ~ 2 eV, which is the largest allowed value from tritium 
beta decay data. The low value of the common mass implied by the WMAP bound leads 
to negligible difference between A 21 and A 21 . This can however be compensated for by 
considering a slightly lower scale for the flavour blind CPT violating mass terms rather than 
the usual Planck scale. 

Acknowledgement: We would like to thank Francesco Vissani for a critical reading of the 



7 




30 
20 
10 

-10 
-20 



f1 



20 40 60 



100 120 140 160 180 



f3 



FIG. 3: P as a function of the phases /i and 
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